General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 






avaiii)bl6 undtr HAS^ spon^rshlp 
interest of esriy •nd wide dis* 
ition cf Earth Resources Survey 

i- ii 'ii ••on 30(1 williout liability 


INVESTIGATION OF REMOTE SENSING TECHNIQUES AS 

B INPUTS TO OPERATIONAL RESOURCE 

MAI^AGEMENT MODELS 


Dr. Fr«d A. Schmwr 
Mr. Robert E. bakson 
RmimI* Sensing ImtHuto 
»$outh Dakfiita Stat* University 
BrookingSf South Dakota^T^S^ 


1976 

Interim Typo n Report far Foried 1 1 Mar 76— 10 Jun 76 





I 


I 


. I ) h I 


TECHNICAL REPORT STANDARD TITLE PAGE 


1. Report No. 2. Government Accession No. 

5 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

5. Report Date 

Investigation of Remote Sensing Techniques as 

July 1976 

Inputs to Operational Resource Management Models. 

6, Performing Organization Cod,e 

7. Author(s) 

8. Performing Organization Report No. 

F.A. SyihjTier, R.E. Isakson 

SDSU RSI-76 05 

9. Performing Organization Name and Address 

10. V/ork Unit No. 

Remote Sensing Institute 


South Dakota State University 

11. Contract or Grant No. 

Brookings, South Dakota: 57006 

NAS5-20982 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Type II Interim 

NASA 

March 11 - June 10, 1976 

Goddard Space Flight Center 


Fred Gordon - Technical Monitor 

14. Sponsoring Agency Code 

Greenbelt, MD 20771 



15. Supplementary Notes 


Participating South Dakota State Government Agencies include the State 
Planning Bureau, the Department of Game, Fish and Parks, the Department 
of Natural Resources Development and the Black Hills Conservancy Subdistrict. 

16. Abstract LANDSAT and aircraft remotely-sensed data are being evaluated by South 
Dakota Government Agencies as potential data sources for present operational 
state programs. Emphasis is placed on procedures which can be easily initiated 
using present state government restraints. The study site is the 824,000 

j hectare (2,035,000 acre) Belle Fourche River Basin in western South Dakota. 

j Evaluation of a NASA-developed LANDSAT CCT rectification program is reported. 

Land use studies center on determining optimum date for Level 1 visual LANDSAT 
interpretations. LANDSAT surface water i nventory capabilities (visual and 
digital) are also being evaluated. LANDSAT CCT's, in conjunction with DIRS have 
been investigated for aspen mapping and crop identification. Drainage analysis, 
using combinations of 1:24,000 USGS Quadrangle maps, high altitude photography 
and LANDSAT imagery, has been conducted and has resulted in an operational 
mapping program for western South Dakota. Entry of resource management data 
into a computerized data base is also under investigation. 


17. Key Words Suggested by Author 18. Distribution Statement 

LANDSAT 

remote sensing 
State of South Dakota 
operational 

data management 

19. Sec urity ClassiL (of this report) |20. Security Classif. (of this page) 

u 




u 


21 .. No.- of Pages 

35 


22. Price 












PREFACE 


This is the fifth quarterly report on an 18-month resource man- 
agement study funded through the NASA LANDSAT foll.pw-on activity.- 
The primary objective of the study is to assist participating state 
agencies evaluate remote sensing (with emphasis on unmanned resource 
satellite data) as a cost-effective data source in operational programs. 
Participating South Dakota agencies include: The Black Hills Conservancy 

Sub-district, the State Planning Bureau, the Game, Fish and Parks De- 
partment, and the Department of Natural Resources Development. With 
the exception of the Subdistrict, each agency has committed half a 
man-year to the project. In this way remote sensing training and 
evaluation of useful data products is enhanced. With NASA high altitude , 
aircraft imagery, low-altitude aircraft data, and ground truth information, 
LANDSAT satellite data are being evaluated as a data source for such 
programs as land use mapping, surface water inventories, aspen mapping, 
and crop identification. Both visual and computerized interpretive 
techniques are being investigated. Also under investigation is a 

NAS'^A-dev'dsloped Digital Image Rectification System (DIRS) which rectifies 

I 

LANDSAT multi -spectral scanner digital tape data using Universal Transverse 
Mercator coordinates as the base grid. The DIRS package will be eval- 
uated for cost effectiveness as well as for the potential of improved 
digital data analyses. Resource data handling is being addressed via 
MAPCLASS, a program which provides a method of data storage and handling 
with output via line printer, color television monitor, and/or CALCOMP 
plotter. Basic project results will provide recommendations on the 
operational use of LANDSAT data' for South Dakota State government. 
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INTRODUCTION 


This is the fifth quarterly report of LANDSAT Follow-On Project, 

NASA Contract NAS5-20982. LANDSAT and aircraft remotely-sensed data 
are being evaluated by three South Dakota governmental agencies: the 

Game, Fish and Parks Department (GF&P), the State Planning Bureau (SPB), 
and the Department of Natural Resources Development (DNRD). Remote 
sensing data and associated interpretive techniques are being studied 
to determine input applicability to present operational state programs. 
Emphasis is placed on procedures which can be easily initiated using 
present state government financial and personnel restraints. The study 
site is the 824,000 hectare (2,035,000 acre) Belle Fourche River Basin 
in western South Dakota. Interpretations reported herein center around 
land use evaluations, surface water mapping and aspen delineation using 
LANDSAT imagery. An MSS CCT rectification program is also being evaluated. 

PROBLEMS 

Late delivery of initial aircraft and LANDSAT data products did 
not allow for tasks to be accomplished according to the proposal 
schedule. It is felt that substantial contributi ons remain and a 
request has been made through the NASA Technical Monitor to extend the 
project. The extension is necessary to ensure completion of the project 
as pr-oposed. Tentative verbal confirmation of the extension has been 
received from the Technical Monitor. 
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ACCOMPLISHMENTS 

Digital Ima ge Rectifica tion System 

One of the major advantages of rectifying LANDSAT CCT data is the 
resultant potential for sequential computer data assessment. The digital 
analysis of rectified CCT's might allow automatic comparison of surface 
water inventories at different time periods. Monitoring of crop re- 
flectance and knowledge of area crops might assist in accurate areal 
crop identification. The potential for automatic mapping and monitoring 
of resources (such as aspen) is well within the possible uses of 
rectified LANDSAT CCT data. Necessary for such work is a method of 
reducing the distortions inherent in the CCT's as supplied by EROS. 

NASA hcS developed a rectification pronram package specifically for 
"correcting" LANDSAT data. The Remote Sensing Institute obtained the 
NASA-developed programming package called Digital Image Rectification 
System (DIRS). DIRS was designed to operate on digital LANDSAT f^S 
imagery and to cartographical ly correct the data to a Universal Trans- 
verse Mercator (UTM) grid system (1). 

Preliminary work was begun during the last quarter (2). Selection 
of approximately 20 ground control points (GCP) per scene is recommended 
(1). The GCP's are geographic or physiographic landmarks which can 
be found on the LANDSAT image and on 1:24,000 USGS Quadrangle maps. 

After the GCP's are selected and the approximate pixel coordinate 
for each GCP is established, DIRS produces a shade print of the GCP 
area. The print is used for verification of the GCP location. The 
pixel '"CO '■u mates obtained from the shade print are used as input 
to the next DIRS step which provides for expansion of the MSS data 
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at the GCP site. These expanded maps are used in conjunction with 
1:24,000 Quad maps to correlate the IJSGS and DIRS GCP coordinates. 

In some areas of western South Dakota there are no 1:24,000 maps avail- 
able. In other cases, the 1:24,000 maps have no UTM coordinates. 

These facts must be kept in mind while selecting GCP's. 

With the coordinate input, DIRS generates data that are used for 
the rectification of an entire scene or portions thereof. To date 
LANDSAT scenes E-21 86-1 7004 and E-222 2-16595 have been processed 
using DIRS. Preliminary evaluation of the recti fi cation/overl ay ac- 
curacies of the program has begun. Initial work with DIRS revealed 
that considerable technician time is required in the processing of 
GCP's. Selection of points, working with the various printouts and 
detemiining coordinates was not a minor consideration for the first 
CCT to be rectified. Subsequent rectification of scenes over the 
same area allowed duplication of GCP's and associated UTM coordinates. 
Based on these procedures approximately 5 man-days are required for 
10 GCP's. As techniques are developed and refined this time estimate 
nay be revised. Based on technician wage of $10/hr the labor costs 
for a scene would be about $400. Computer rectification costs for 
CCT E-2222-16595 are sui.snarized in Table 1. The "Definition of Rec- 
tification parameters" step calculates the data necessary to actually 

rectify specific areas within a scene or the entire scene. In an 

2 2 

aspen mapping study, a 93 km (36 mi ) area was rectified and line 

2 2 

printed for an average cost of $. 07/km ($.17/ini ). 

Total costs can be summarized as in Table 2. One of the major 
advantages of rectifying LANDSAT data is the resultant potential for 
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sequential computer data assessment. Necessary for such work is an 
overlay capability of sequential LANDSAT data. Such overlay capability 
will allow monitoring of specific points on the ground, on a temporal 
basis, by permitting pi xel-for-pi xel overlay in a time-sequential 
period limited by cloud conditions and LANDSAT orbit configurations. 

To evaluate DIRS output for such requirements, a series of geographic 
accuracy test points have been selected. For inclusion as a test 
point, an object must be identifiable both on the 1:24,000 Quadrangle 
map and on the DIRS-rectified output. A point is selected and the 
latitude and longitude is determined from the base map. A four minute 
by four minute window surrounding the point is extracted from the 
CCT and processed using DIRS. Knowing the window boundary latitude 
and longitude and spotting the test point on the printout allows 
calculation of the DIRS-generated latitude and longitude by simple 
scale measurements and ratio! ng. The geometric accuracy test points 
thus far evaluated are listed in Table 3. Additional test points 
are planned. To assist in evaluating the information contained in 
Table 3, the differences in combinations of latitude/lonqi tude are 
presented in Table 4. 

It should be reiterated that identical GCP's and UTM coordinates 
were used for the July and September CCT's. There are two basic reasons 
for this: 1) decreased costs involved in GCP processing and 2) ability 

to enhance potential overlay capabilities of the rectified output. 

It was felt that the spatial accuracy of DIRS coordinates as compared 
to base map coordinates was not as important, from a sequential CCT 
analysis point of view, as was overlay capabilities of various LANDSAT 


6 



I 


I 



» 


. 

- 

- 1 

j 



0 

UO 

00 1 

>- 



0 

LO 

0 



cn 

“ 

— 

- 

<c 


c 

LT> 

r-^ 


s 


0 

CO 

c>u 

0 



_J 

0 

0 

0 

O 



CO 

CO 

CO 

o 



0 

0 

0 

et 



' 



o 

►—1 

Q- ^ 

, 3 ^ y 


- 

- 

- 

•X 

L /0 


CO 

>x> 

LO 

CL 




CVJ 

CO 

s 


4 -> 

» 

“ 

•• 


03 

00 

CO 

CTi 

C 3 OO 


_j 

r-* 

r— 

0 

O LlI 



0 

0 

0 

LiJ 1 — 





«=t* 

C 3 eX 




<d“ 

Cj“ 







Q 

LlJ O 



- 


- 

t— ci: 



CNJ 

CO 

cn 

g§ 



0 

LO 

0 


cn 

— 

“ 

“* 

LU O 


c 

LD 

r— " 

1 



0 

CO 

C^U 

0 

LU Ou 


-J 

0 

0 

• 

C 3 ^ 



00 

0 

CO 

0 

CO 1 

0 




f— 

r— 








CO Ln 

'3 ^ 





o; cc 




— 


♦ CQ 





1 

Q 




CNJ 

CO 

o 


+J 

— 

- 

— 1 

U_ 1 — 


03 

00 

CO 

S ! 

o 


_J 

r— 

r— 

<=> 1 

Q 



0 

0 

0 

LU LU 






o a: 



1 ^ 



K- Q. 
— < 2 : 



- 

r 

“ 

o o 



0 

LO 

CNJ 

Z CJ 



0 

0 

r— 

o 


01 

— 

— 

• 

1 lo 


c 

LO 

CNJ 


ct 


0 

CO 

cvj 

0 

Q 


u 

0 

0 

c 

^ OO 

CL 


CO 

CO 

CO 

< 1 — 

fT 3 


0 

0 

0 

z 

CO s: 


r— 

r— ' 



0 





Q O 

CO <u 





Z 3 Q- 

n) CO 


— 

— 

— 

h- 

ro 


CNJ 

VO 


I 1 1 

CO 



CM 

CO 

H- OO 


+J 

— 

— 

“ 

cc UJ 


03 

00 

CO 

cn 

—I 1 — 


—1 

r— 

r— 

0 




0 

0 

0 







CO 

LU 






0 





_j 






CO 

-*-> 





C 3 : 

LO U 


r— 

CM 

CO 

h- 

O) c 
1 — 






0 






cu 






7 


I 


I 


I 


I 


I 


Q 

UJ 


O 

LU 

q; 


oo 

ai 

»— I 

Q 


Q 

z: 

<c 


Q. 

c 


UJ 

(/I 

<c 

CQ 

q: 

o 


< 

(— 

c 

Q 


UJ 

t— 


a 

Cd 

o 

o • 

O 


UJ I— 

> o 
1—1 o 


r 

§ 


*a: 
C CO 
Q. Q 


o < 

O _J 


UJ 

ja 

<c 



8 


f 


scenes. Use of identical GCP's and associated UTM coordinates should 
~a/imize the overlay efficiency of DIRS. 

Again, Table 4 is not complete as a larger number of test points 
regain to be examined. The averages are not of much statistical validity 
considering the small sample, but they may indicate a trend in accura- 
cies. Considering the 70 m x 70 m DIRS pixel dimension, the average 
errors in overlaying July and September data appear encouraging (40 m 
latitude and 29 m longitude) or, less than the specified DIRS pixel. 
Continued analysis of other CCT’s, additional GCP's, geometric accuracy 
test points and costs should allow basic conclusions to be drawn re- 
garding DIRS applicability. 

DIRS output is not north oriented, i.e. line printer and plotter 
maps are rectified but remain rotated approximately 12 degrees west 
of north. This is not a serious problem in classification work but 
it is a restraint to the potential automatic entry of data into a 
data base like MAPCLASS. Development of a rotation routine is desirable 
but not within the scope of this project. 

Land Use 

All participating state agencies have expressed interest in eval- 
uating various land use data sources. Game Fish and Parks (GF&P) and 
the Department of Natural Resource Development (DNRD) are under leg- 
islative mandate to provide resource inventories of each of the state's 
16 river basins every four years. One of the most important resources 
in South Dakota is the land itself. Detailed information regarding 
its use and the spatial distribution of the various uses is considered 
as significant irjut for these four-year surveys. The State Planning 


9 



I 


Bureau (SPB) is preparing a landuse map of the state using digital 
analysis of LANDSAT MSS CCT data. The Bureau is interested in eval- 
uating the digital analyses v/ith results obtained from other data 
sources . 

Historically, land use information has been obtained from a va- 
riety of sources including the Conservation Needs Inventory and Crop 
and Livestock Reporting Service Reports. Usually the data are not 
subdivided into river basins but into political subdivisions such as 
counties. This procedure results in estimating errors as the data 
are manipulated to a river basin status. 

To assist in the evaluation of the accuracies of LAMDSAT inter- 
pretations, an area which contained a good selection of cloud-free 
LANDSAT and high altitude imagery was chosen. The site, a ten-township 
230,400 acre (93,279 hectare) area, is located north and east of Sturgis 
in the center of the basin. 

Two dates of high altitude coverage are available for the site. 
Interpretations of color infrared transparencies produced the agri- 
cultural percentages listed in Table 5. The discrepenci es between 
the two interpretations will be discussed later. 

TABLE 5. HIGH ALIIFUDE AIRCRAFT INTERPRETATION 


Data Source 

Date 

Scale Interpreted 

Percent Agricultural 
Land 

Hi gh-al ti tude 




Color IR 

25 Jun 75 

114,506 

27 

High-al ti tude 




Color IR 

23 Sep 75 

125,880 

22 
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A variety of LANDSAT data are available and i nterpretations have 
been conducted as listed in Table 6. Sufficient data were available 
to allow representation of each of the months May through November. 
Additional input from DNRD and GF&P interpreters is expected. 

Optimum interpretations would result in 100 percent of the RB-57 
agricultural lands (in terms of total area) being interpreted. However, 
selection of optimum dates for interpretation should not be based 
on such information alone. Spatial accuracies were determined for 
the interpretations as indicated and include the optimum (gross accuracy) 
i nterpretation available for May through November. 

The spatial accuracy data are listed in Table 7. Rangeland is 
consistently interpreted as rangeland approximately 90 percent of 
the time throughout the season. Agricultural land i nterpretations 
are much more variable. It appears from the data in Table 7 that 
a prime date for interpreting agricultural land from LANDSAT for western 
South Dakota, is in the mid June to late July growing season. 

Referring to Figure 1, a Belle Fourche Basin crop calendar, the 
mid June to July period intersects most crops in LANDSAT-recognizabl e 
stages of development. Small grains are heading and ripening, iden- 
tifiable as red or yellow patterns, respectively on a color composite. 
Corn and sorghum can be identified as either a bare field shortly 
after planting or as a bright red color as the plants grow and mature. 
Alfalfa is identifiable with a very bright red signature. Wild hay 
areas, which, for this interpretation would be considered agricul tural , 
ar(' interpreted only after cutting and then as a pattern within the 
predomi nantly range areas. 
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TABLL 6. COLOR COMPOSITE INTERPRETATIONS FOR lO-TOWNSHIP AREA. 


Date 

scale 

Source* 

Ag. land as a 
percent of interpreted 
RB-57 Ag. land 

24 May 

1:125,000 

E-P (DNRD) 

65 

24 May 

1:125,000 

E-P 

104** 

24 May 

1:125,000 

R-P 

123 

12 June 

1:125,000 

E-P 

93** 

18 July 

1:125,000 

E-P 

134 

27 July 

1 : 250, 000 

ZTS 

126 

27 July 

1:125,000 

ZTS 

130** 

27 July 

1 : 250, 000 


130 

27 July 

1:250,000 

E-P (DNRD) 

56 

27 July 

1:125,000 

E-P 

111** 

27 July 

1:125,000 

R-P 

102** 

27 July 

1:125,000 

E-P (DNRD) 

TOO** 

5 August 

1 :125,000 

E-P 

95** 

23 August 

1 :125,000 

E-P 

96** 

1 

1 September 

1:125,000 

E-P 

100** 

1 September 

1 :125,000 

R-P 


October 

1:125,000 

E-P 

51 

November 

1:125,000 

E-P 

73 


* E = EROS Color Composite 
R = RSI Color Composite 
P = Print 

ZTS = Zoom Transfer Srope 

** Accuracy %'s calculated 
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Figure 1. Crop calendar for the Belle Fourche River Basin with optimum LANDSAT Level 
Land use window indicated. 
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Interpr'ftation accuracies for three sets of black and white inter- 
pretations {> 1 ) are presented in Table 8. Additional interpretations 
will be conducted prior to project completion. 

Comparisot! of June RB-57 and September U-2 interpretations resulted 
in an interesting observation. If it is accepted that the two initial 
data sets were of comparable quality, e loss of five percent in the 
overall agricul tural areas is noticed in the September interpretation. 
Spatial comparison with June RB-57 data is seen in Table 9. It is 

TABLE 9. spatial COMPARISON OF SEPTEMBER HIGH ALTITUDE AIRCRAFT DATA 
WITH JUNE HIGH ALTITUDE AIRCRAFT DATA 


Ag as Ag 

Ag as Range 

Range as Range 

j 

Range as Ag* 

76.4 

23.6 

■ 97.1 ■ ! 

! ■■■ ;:9 


* Includes 720 acres of new agricultural land as rangeland 
interpreted as agricul tural . 


suggested that the 24 percent agricultural land classified as rangeland 
from the September interpretation is a result of land cover changes 
as the season progresses. There are a number of interpretive guides 
used to assist in classifying land into the agricultural category 
from high altitude imagery. Color, texture and patterns all assist 
and if classification remains in doubt, the apparent resolution of 
the film allows recognition of detailed information such as haystacks, 
cutting and cultivating patterns, ecc. In the June interpretation , 
all of these interpretive devices could be used. By September, color 
differences are reduced substantially and the ititerpretation was heavily 
based on other indicators. Certain pasture areas in particular were 
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unidentifiable without the color differentiation present in the June 
imagery. As LANDSAT interpretations must rely soley on the color/ pattern 
signatures, it is reasonable to expect an increase in LANDSAT-based 
classification of agricultural land as range (Table 7) as the season 
progresses. It is important, then, to obtain ground truth (i.e. aerial 
photography) during the time frame also designated for optimum satellite 
i nterpretation. 

Surface Water 

Analysis of interpretive scales and data sources was initiated 
to provide a basis for recoimiendations regarding LANDSAT as a surface 
water inventory tool. While surface water is a category in Level I 
land use, it is handled separately here. To be classified in a Level I 
land use classification, the surface area must equal or exceed 40 acres (16 
(16 ha). In western South Dakota there are numerous bodies of water 
which are considerably less than 16 hectares. The number and distrib- 
ution of these smaller bodies of water are important in evaluating 
the general water (and soil moisture) conditions throughout the basin. 

The small water bodies are also important to western South Dakota 
fislieries management, wildlife and livestock production. To ifiventory 
only 16 ha and greater surface water bodies, as specified in Level 1 
Land Use, would place serious restrictions on the reliability and 
usefulness of the surface water inventory. For these reasons, surface 
water is separated from the Level I Land Use category and the surface 
water inventory capabilities of LAMUSAT are separately investigated 
in detail . 


17 


t 


I 


1 


The same 10-township area used in land use analysis is being 
used for surface water studies and can be seen in Figure 2 . Initial 
interpretations were concerned with a strictly numerical accounting 
of bodies of water within the 10-township area. A variety of dat.^ 
sources were employed and are 1 isted in Table 10. A number of inter- 
pretations have been added since the last guarterly report (2). 

An ASCS photo interpretation was prepared by DNRD. This source 
was included because DNRD has interpreted surface water from such 
photographs for the entire state and the Department is interested 
in comparing the results of their present procedures with what might 
be available via LANDSAT. The comparative inventory value (4) of 
the ASCS interpretation is better than most LANDSAT values. But, when 
the interpretation is more closely analyzed (within the smaller area 
of Fig. 2), the apparent advantage over LANDSAT is not so pronounced. 

Referring to Table 11, the ASCS i nterpretation resulted in 
77 percent identification of water bodies over 1.1 acres (.45 ha). 
Interpretations of f4SS7 data produced "above 77 percent" figures 
in all but four LANDSAT Visual i nterpretations . The ASCS data 
were more adequate than LANDSAT for interpreting the less than 
1.1 acre water bodies. The apparent resolution of the imagery is 
offered as explanation. The des'^repencies between LANDSAT and 
ASCS i nterpreta t ions of > .45 ha data might be the result of 
interpretive variabi 1 i ti es , the fact that the ASCS data are nearly 
ten years old, and, the fact that the ASCS imagery was collected 
in vluly, a potentially dry period. 
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f-'qure 2. Surface water interpretation of the Relle Fourche River 
Basin as interpreted by OrjRD personnel from snrinn 1975 
LANOSAT MSS 7 imaqery. Visual print i nterpretation was 
conducted at 1:250,000. The larner rectanoular area is a 
ten-township site interpreted for water bodies from a 
variety of data sources. The smaller outlined area was 
interpreted with more detailed evaluation of LANDSAT 
inventory capabilities. 




TABLE 10. SURFACE WATER INTERPRETATIONS FOR lO-TOWNSHIP TEST AREA (Townships 6 and 7 North) 
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A CCT of scene 2186-17004 was rectified and level -si iced for 
water using SPB-developed programming packages. Window selection 
allowed processing of T6N only of the ten township area depicted in 
Figure 2. Color infrared and LANDSAT print interpretations were ob- 
tained for T6N and are listed along with the digital data in Table 10. 

The visual and digital procedures appear quite similar in relation 
to the base hi gh -al ti tude data. Detailed analysis of the smaller 
area of Figure 2 (Table 11) indicates threshold values were not low 
enough to allow detection of any less than .45 hectare water bodies. 
Program capabilities are big developed which would allow detection of 
these smaller bodies . 

Within the smaller area of Figure 2, 49 bodies of water interpreted 
both on RB-57 ground truth and LANDSAT MSS7 were used in estii ation 
of water surface area fr'om visual LANDSAT interpretations. The resultant 
data are found in Table 12. Area measurements of high altitude aircraft 
imagery were accomplished using 1:24,000 enlargements and a grid method. 
LANDSAT i nterpretation areas were estimated by "■’‘sual comparison of 
the interpreted dot with the 40 acre grid cell used in digitizing 
Tor MAPCLASS. From the data it appears that surface water area is 
overestimated by .8 - 1.2 hectares per average water body using visiial 
LANDSAT interpretation methods. 

Table 13 provides low altitude aircra<"t data which is indicative 
of the seasonal variability of surface water referred to in the lANDSAT 
interpretations of the last quarterly report (2). For the interpreted 
sections of Table 13, both the total number of water bodies as well 
as total surface ar’ca has decreased from June to Sept niber. These 
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data were further refined by subdividing the surface area measurements 
into categories as listed in Table 14. For the sections under con- 
sideration, the largest number of water bodies are less than .2 ha 
with a majority being less than .5 ha. The importance of inventorying 
these smaller water bodies and the ability of LANDSAT to conduct such 
an inventory will be addressed during the course of this project. 

? i^op Id enti f ic ation 

LANDSAT color composite data o*" test area number six (3) were 
printed to a scale of 1:125,000. An attempt was made to identify 
fields within the ground truth area. Considerable visual inspection 
allowed the location of only the largest of the area fields. Within 
the Belle Fourche Irrigation District (which is where ground truth 
was collected) fields are not large, average sizes of 4.5 ha are not 
unconinon. These size fields were not easily recognized by visual 
i nterpretation of 1:125,000 color composite LANDSAT data. 

An attempt at digital analysis of sequential CCT's using DIRS- 
rectified LANDSAT data again resulted in an inability to identify 
known fields. Shade prints of MSS5 and MSS7 for 27 July and 1 September 
were used in the analysis. Boundary-detection techniques co'jld be 
employed to assist in field locations. However, from an operational 
point of view, it appears accurate I.ANDSAT-based crop identification 
on a river basin basis is not feasible at this time and additional 
work will not be emphasized. 
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TOTALS 39 19.4 
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Dr ainage Systems 

An interest was expressed in drainage networks and basin/sub-basin 
delineation. An accurate d.ainage map of the Basin is not in existence, 
so an effort was pi t forth to demonstrate and develop a procedure 
to produce such a map. The basic premise was to work from available 
1:24,000 USGS Quads. The drainage information on quadrangle maps 
was considered as the optimum readily available data. Where quads 
were not available, high altitude aerial photography was designated 
as prime data source. Where neither quad nor aerial photography was 
available, LANDSAT imagery was prime data source. 

The developed procedure is as follows. Transfer drainage infor- 
mation from 1:24,000 Quadrangle maps to velum. Where maps are un- 
available, a ZTS and high altitude color infrared transparency are 
used to develop requisite drainage data. In areas .-.here no aerial 
photography or drainage maps are available, drainage informiation is 
obtained from a combination of interpretation of winter black and 
white and summer color composite sccfies. I.ANDSAT winter scenes ac- 
centuate certain drainage features via the combination effects of 
snowcover, and sunangle. In very flat areas the winter data is difficult 
to interpret. For these areas, the color composite provides drainage 
analysis indicators such as vegetation along stream and soil erosion 
signatures. A combination of the two seasons allows for the generation 
of a LANDSAT-based drainage map in areas devoid of other data. 

All interpicted die inage data are photographed and printed to 
the same scale for preparation of a controlled mosaic. Tfie mosaicked 
dr-iinagc map is then drafted as the final drainage map. 
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Continued analysis will allow comparison of LANDSAT-devel oped 
drainage data with that of 1:24,000 Quadrangle maps. Temporal overlays 

r> 

using (Model 6040 PT) is also contemplated. 

The methods described above have resulted in similar work being 
conducted by RSI for the Sixth Planning District for seven western 
South Dakota counties. In addition to basic drainage mapping, basin/sub- 
basin delineations are interpreted and surface water is interpreted 
from LANDSAT MSS7 imagery. The resultant products, county-wide hydrology 
maps, are being prepared in this operational situation as a direct 
result of LANDSAT Follow-On Project developments. 

Fl oodplain A na 1^ sj s_ a n^ S t r ij) De ve 1 ppnejT^ 

Planning Bureau personnel have historically developed floodplain 
maps from existing contour maps. Using this procedure, deve1opi;ent 
of flood high water mark has been an intuative procedure based on 
contour lines. The Bureau is in the process of evaluating high altitude 
aerial photography for flood plain analysis. A region surrounding 
Belle Fourche, South Dakota is being mapped For floodplain using two 
data sources: 1:24,000 USGS Quadrangle Maps and high altitude imagery. 

Extensive private development along the "strip" area between 
Sturgis and Rapid City has prompted the State Planning Bureau to map 
the entire area for land jse. A map of the Blackhawk Quadrangle is 
seen in Figure 3. The information was interpreted from 1975 high 
altitude transparenc ies using a Zoom Transfer Scope. Low altitude 
RSI imagery was viewed stereoscopically to veriFy any questionable 
classifications. An interpretation of 1974 imagery was available 
via NASA Office of University Affairs, Grant No. NGL 42-003-007 and 
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Fifjure 3. I afuJuso dlotuj St.urcji s-Rupi d City Strip area as interpreted 
ffom 197S aircraft data. 





is seen in Figure 4. Note the extensive development in the southern 
portion of the 1975 map as compared to the 1974. Development is so 
rapid that it is becoming a problem for local government officials. 

Remote sensing offers a potential method for nwnitoring the development. 

As pen 

An aspen resource inventory of the Black Hills area has assumed 
added importance to GF&P personnel. The Department was involved in 
a study which indicated pulverized aspen to be a good cattle feed 
supplement. Black Hills aspen groves may now become an asset. Potential 
development includes planned harvesting of aspen for use as a feed 
supplement. Development of a method for locating and monitoring aspen 
groves resulted in GF&P initiation of this aspect of the project. 

An area of representati ve forest types was located using high 
altitude aircraft imagery. Visual interpretation of July LANDSAT 
color composite data was deemed impractical (4). Digital analysis 
of CCT data proved mor'e encouraging (2). Subsequent digital analysis 
using RSI's K-class package in conjunction with NASA's DIRS program 
and a CALCOMP plotter have allowed geometrically corrected LANDSAT- 
based aspen maps. The "aspen as aspen" classification is at least 
50 percent better using LANDSAT data than is available from existing 
US Forest Service Maps of test township T4NR5F. In certain areas, 

LANDSAT digital analysis identified aspen groves which were not indicated 
on the high altitude i nterpretation. 
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Me pci a ss 

Considerable effort has been expended iti the evaluation of entering 
remotely-sensed resource i nfonTiat''on into a computer-based data system. 
MAPCLASS, as the program is called, allows for manual coding of change- 
points rather than coding every cell in a grid system. The program 
also allows the computer to "overlay" and manipulate numerous combi- 
nations of data entered into the data base. Output via a Calcomp 
plotter permits maps to be generated at any desired scale. 

The basic cell is 16.2 hectares (40 acres). This size was selected 

because it is small enough to obtain desired detail, but large enough 

to keep digitizing times reasonable. The grid system ilself could 

have been based on a variety of coordinate systei'is. Two systems cJ^Jpeared 

to be the most prnctical : Latitude/Longitude and Range/ Township. 

In selecting one of these two systems as the grid base, it was thought 

the Latitude/Longitude grid system would be easier to establish; however 

few South Dakotans are familiar with Latitude/Longitude as a ground 

reference system. South Dakota is subdivided into townships by a 

Range and Townsfiip gr id system and because the nomenclature is widely 

used and understood, the decision was made to evaluate a data-base 

grid system t)ased on Range and Township lines. A 1:250,000 base map 

of the Belle Fourc he basin was prepared from 1:250,000 USGS Quadrangle 

maps. Based on Range correction lines and adjusted township 1 ’ les , 

the basin was divided into eigtit sections varying in size frjit; eight 

2 2 

to 30 townships (03.3 km , or 36 mi per township) each. These eight 
sections ar'r.* essentially square in themselves and are each the basis 
of an individual grid systrnn. The digitizing grid system was generated 
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jsing the Ca'tcomp plotter and asso( iati.'d software. With the grid 
system established, it remaitiod to ('Vcihiatt* the digitizing procedure 
and the practicality of the results. 

A variety of data have been entered into the MAPCLASS data base. 
Coded information include: soil association, slope, land use, sub- 

basins, Black Hills National Forest Boundary, counties, Range-Township 
lines, and surface water. Demonstration overlays of the various data 
have been prepared for evaluation by state agencies. Cost information, 
based on coding and processing the above data sets, is being prepared 
for inclusion in a cost effective analysis by partici Dating agencies. 

FUNDS EXPENDED 

Total funds expended through May 31, 1976: $86,440.14. This 

does not include costs incurred by state agency participants as they 
invoice on a quarterly basis. 

DATA USE (as of May 28, 1976): 

Value of Data Allowed - $15,144.00 
Value of Data Ordered - $11,316.00 
Value of Data Received - $9,368.00 

AIRCRAFT DATA 

High and low altitude aircraft imagery have o.^en used ext(nsively 
as a ground truth supplement for surface water, land use, and aspen 
management studies. The aircraft data have been used as primary remote 
sensing data sources for' floodplain mapping and level II landuse 
'tipping. Cloud conditions over certain ground-truthed areas, pr-eviously 
designated as intensive study sites (3), has necessitated r-elian<e on 
high altitude aircraft iV-agery as a j.-ta hise. 

.J 


I 


I 


I 


I 


I 


PROGRAM FOR NEXT REPORTING PERIOD 

Further aspen analysis will quantify mapping accuracies. Additional 
land use interpretations using black and white LANDSAT imagery is planned. 
Development of input for DNRD's Belle Fourche River Basin Inventory is 
expected. Continued evaluation of DIRS is planned to continue inves- 
tigation of relative accuracies involved. Cost effectiveness will be 
investigated. Input from State agencies is expected and will provide 
participi tory information as well as recommendations and conclusions 
on the use of remote sensing data within each agency. 

CONCLUSIONS 

Continued LANDSAT analysis by and with participating state agencies 
has demonstrated to them the relative ease with which certain infor- 
mation can be obtained from the imagery. Comparisons with other 
available data, including the previous data source for a state-wide 
inventory, demonstrate the value of LANDSAT data as it relates to an 
accurate accounting of western South Dakota surface water. Mid June- 
July LANDSAT data appears an optimum source for Level 1 Landuse inter- 
pretations. Digital LANDSAT CCT analysis for aspen has produced 
maps 50 percent more accurate tlian any present typing maps of the 
test area. LANDSAT data have been successfully used in a drainage 
network mapping procedure which has resulted in an operational program 
for seven western South Dakota counties. Relatively simple and 
straightforward procedures have been employed in all analyses and the 
results demonstrate the potential value of the data. 
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